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Abstract 



We examine a result of Basor and Ehrhardt concerning Hankel and Toeplitz plus 
Hankel matrices, within the context of the Riordan group of lower-triangular matrices. 
This allows us to determine the LDU decomposition of certain symmetric Toeplitz plus 
Hankel matrices. We also determine the generating functions and Hankel transforms 
of associated sequences. 
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1 Introduction 

In [1] Basor and Ehrhardt studied the transformation 



defined for sequences {cin}'^=-oo the context of relating the determinants of certain Toeplitz 
plus Hankel matrices to the determinants of related Hankel matrices. 

In this note, we shall study an equivalent transformation, which we will construct with 
the aid of Riordan arrays [9]. We call this the B-transform. We shall then use our results to 
examine the LDU decomposition of the resulting Toeplitz plus Hankel matrices. 

In the next section, we shall detail the notations that will be used in this note, and 
give a basic introduction to the relevant theory of Riordan arrays. We shall follow this 
with a section which defines the B-transform, studies some of its properties, and shows its 
equivalence the Basor and Ehrhardt transform. In particular, we derive an expression for 
the generating function of the image sequence, which for instance allows us to determine 
the Hankel transform of image sequence in many cases. A final section then looks at the 
LDU decomposition of the related Toeplitz plus Hankel matrices, with examples involving 
Riordan arrays. 
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2 Notation and basic Riordan array theory 



Although many of our results will be valid for sequences a„ with values in C, we shall in 
the sequel assume that the sequences we deal with are integer sequences, a.„ G Z. For an 
integer sequence a^, that is, an element of Z^, the power series f{x) = Yl^=o'^kX^ called 
the ordinary generating function or g.f. of the sequence. a„ is thus the coefficient of in 
this series. We denote this by a„ = [x'^]f{x). For instance, F„ = [x"] ]^_^_^2 is the n-th 

Fibonacci number, while C„ = [x"] ""^"^^"^^ is the n-th. Catalan number. We use the notation 
Qn _ I'^n.j]^ fQ]. sequence 1, 0, 0, 0, . . . Thus 0" = [n = 0] = 5nfl = (°). Here, we have used 
the Iverson bracket notation [3], defined by \P] = 1 if the proposition V is true, and [P] = 
if V is false. 

For a power series f{x) = ^^q*^"^" with /(O) = we define the reversion or composi- 
tional inverse of / to be the power series f{x) such that f{f{x)) = x. We sometimes write 
/ = Rev/. 

The Hankel transform [6] of a sequence a„ is the sequence hn = \ai+j\^j^Q. If the sequence 
an has a g.f. that has a continued fraction expansion of the form 

Oo 



1 — aox 



P2X 

1 — aix 



1 — a2X — 



1 

then the Hankel transform of a„ is given by [5] 



n 



k=l 



The LDU decomposition of Hankel matrices has been studied in [7, 8]. 

Some of the lower-triangular matrices that we shall meet will be coefficient arrays of 
families of orthogonal polynomials. General references for orthogonal polynomials include 
[2, 4, 11]. 

y4* will denote the transpose of the matrix A, and we will on occasion use A - B to denote 
the product of the matrices AB, where this makes reading the text easier. This also conforms 
with the use of "■" for the product in the Riordan group (see below). 

The Riordan group [9, 10], is a set of infinite lower-triangular integer matrices, where 
each matrix is defined by a pair of generating functions g{x) = 1 -|- gix + g2x'^ + ■ ■ ■ and 
f{x) = fix + /2X^ + ■ ■ ■ where /i 7^ [10]. We assume in addition that /i = 1 in what 
follows. The associated matrix is the matrix whose i-th column is generated by g{x)f{xY 
(the first column being indexed by 0). The matrix corresponding to the pair / is denoted 
by ((?, /) or lZ{g, /). The group law is then given by 

(^7, /)-(/i,/) = (^?, /)(/!, /) = (^?(/.o/),/o/). 

The identity for this law is / = (l,x) and the inverse of {g^f) is {g,f)~^ = (1/(5' o /),/) 
where / is the compositional inverse of /. 
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A Riordan array of the form {g{x),x), where g{x) is the generating function of the se- 
quence Qn, is called the sequence array of the sequence a„. Its (n, k)-th term is an-k- Such 
arrays are also called Appell arrays as they form the elements of the Appell subgroup. 

If M is the matrix [g, /), and a = (ao, ai, . . .)' is an integer sequence with ordinary generat- 
ing function A (x), then the sequence Ma has [9] ordinary generating function g{x)A{f{x)). 
This result is often called "the Fundamental Theorem of Riordan arrays". The (infinite) 
matrix [g, /) can thus be considered to act on the ring of integer sequences by multi- 
plication, where a sequence is regarded as a (infinite) column vector. We can extend this 
action to the ring of power series Z[[a:]] by 

{g, /) : A{x) ^ {g, /) • A{x) = g{x)A{f{x)). 

Example 1. The so-called binomial matrix B is the element {jz^, yz^) of the Riordan 
group. It has general element (^), and hence as an array coincides with Pascal's triangle. 
More generally, B"* is the element (jt^^, j^^) of ths Riordan group, with general term 
("')m"'~'^. It is easy to show that the inverse B""^ of B"* is given by (-rr — , -rr — )• 



For a sequence ao, 01,02, . . . with g.f. g{x), the "aeration" of the sequence is the sequence 

2x 



Oo, 0, oi, 0, 02, . . . with interpolated zeros. Its g.f. is g{x^). We note that since c{x) — ^ 



has the well-known continued fraction expansion 

c{x) = 



X 



1 



X 

1 



1 



c(x^) has the continued fraction expansion 



c(x^) = . (3) 



x'^ 



x"^ 



The aeration of a (lower-triangular) matrix M with general term rriij is the matrix whose 
general term is given by 

1 + (-1)^-^' 

f^i+i_ i-j ; 
2 ' 2 Z 

where m^j is the i,j-th element of the reversal of M: 

In the case of a Riordan array (or indeed any lower triangular array), the row sums of the 
aeration are equal to the diagonal sums of the reversal of the original matrix. 
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Example 2. The Riordan array (c(x^), xc(x^)) is the aeration of the Riordan array 

(c(x), xc(x)) = (1 — X, x(l — x))~^. 

Here 



c{x) 



1 - y/l-Ax 

2x 



is the g.f. of the Catalan numbers. Indeed, the reversal of {c{x),xc{x)) is the matrix with 
general element 

[k<n+l] 



k 



n + 1 



which begins 



/ 1 

1 1 

1 2 2 

1 3 5 5 

1 4 9 14 14 

1 5 14 28 42 42 



V : : : 

Then (c(x^), xc(x^)) has general element 



and begins 



n + l\ fc + ll + (- l) 
^ jn + l 



/ 1 

1 

10 10 

2 1 

2 3 1 

5 4 1 



n—k 



We have 



v 



X 



1 + x^ 1 + x^ 



3 The B-transform 



We let 



1 — X 



X 



1 + X2 ' 1 + x2 
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where we recall that the Riordan matrix {g, f) is the lower triangular matrix whose k-th 
column has generating function g{x)f{x)^ , for suitable g and /. Then L has (n, k)-i]i term 

n 

^^,Yq-^) is the coefficient array of the generalized Chebyshev polynomials 

1. 



and = [l 
defined by 

Pn{x) = xPn-l{x) - Pn-2{X), Po{x) = 1, Pi(x 

The matrix L begins 



X 



( 1 

















1 


1 














2 


1 


1 











3 


3 


1 


1 








6 


4 


4 


1 


1 





10 


10 


5 


5 


1 


1 



V 

We now define B to be the matrix 



This matrix begins 



L • (1 + x.xf 



( 1 


1 














1 


2 


1 











2 


3 


2 


1 








3 


6 


4 


2 


1 





6 


10 


8 


5 


2 


1 


10 


20 


15 


10 


6 


2 



V ■■ 



{i + x,xy 



xY is 


given 


by 




















/ 1 


1 











. 


.. \ 

















1 


1 








. 






















1 


1 





. 

























1 


1 


. 




























1 


1 . 































1 . 
















v • 












•• / 














/ 1 














. 


.. \ 




/o 


1 

















1 











. 












1 

















1 








. 















1 

















1 





. 




+ 














1 

















1 


. 





















1 

















1 . 

























(4) 
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we see that the (n, k)-th element of B is given by 

Now let {a„}„>o be a sequence. We define the B transform of a„ to be the sequence 
given by 

n+l 
k=0 

where B = (6n,fc)n,fc>o- 

Example 3. The B-transform of the Fibonacci numbers is given by 
This sequence starts 

1, 3, 7, 17, 39, 91, 207, 475, 1075, 2445, 5515, .... 
It has the interesting property that its Hankel transform is (—2)"'. 
Proposition 4. We have 

®=(rT^'rf^) 

where T is the matrix 



Proof. We have 



Now 



T= ( Y^,x] ■{l + x,xy. 



L - {l + x,xy = L - i ■ (y^^' -{l + x, xY. 



L-l^xV = " W ^ x~' 

'l-x' / Vl + a;2'l + x2y \l-x' 

1 \ / 1 — X X 



1 — a;' / yl + x^'l + x^ 

1 ^ -1 

1 X 

1 + ' 1 + 
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We recall that the matrix {j^, jf^) ^ = (c(x^), xc(x^)), where c{x) 
of the Catalan numbers, has general element 



2x 



is the ff.f. 



n + l\k + ll + (-1) 



n—k 



n + l 2 

In addition, (j^^, x) ■ (1 + x, x)* is the matrix T given by 

/ 1 1 ... \ 

12 10 

12 2 10 

1 2 2 2 1 

1 2 2 2 2 1 

1 2 2 2 2 2 



T 



v 



Now note that 



T 



/ 1 














. 






/o 


1 











. 


.. \ 


1 


1 











. 









1 


1 








. 




1 


1 


1 








. 









1 


1 


1 





. 




1 


1 


1 


1 





. 




+ 





1 


1 


1 


1 


. 




1 


1 


1 


1 


1 


. 









1 


1 


1 


1 


1 . 




1 


1 


1 


1 


1 


1 . 









1 


1 


1 


1 


1 . 




v • 
















v • 















and hence the action of T on a sequence a„ is to return the sequence with n-th term equal 
to 



n+l 



^ flfc + ^ ctfc = 2 ^ Ofe + a„+i - 



ao. 



fc=0 



k=l 



k=0 



Thus we have 
Proposition 5. We have 



k=0 



n + 1\ A; + 1 1 + (-1) 



n--k 



k+1 



n + l 2 ^^"^ ■ ^ 
3=0 i=i 



or equivalently, 



/n + 1\ A; + 1 1 + (-1)""^ A 
bn = 2^( ^ {2^a,+ak+i-ao). 



k=0 



j=0 



In the following, we will be interested in determining the g.f. of the image of a^. We 
have the following result. 
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Proposition 6. Let f[x) he the g.f. of an- Then bn = Y12=o ^-n^kf^k has g.f. given by 



1 X \ ^ f {1+ x)f{x) - ao 



1 + + \ x{l — X 

Equivalently, the g.f. of bn is given by 

1 - X X / {l+x)f{x) - ao \ 

1 + x^'l + x^J '[ X J' 

Proof. The result follows from the fact that the generating function of Ylk=o + Ylt=i 
is given by 

1 1 ffix)-ao\ (l + x)/(x)-ao 



1 — X 1 — X \ X J x{l — x) 



Corollary 7. The g.f of bn = El =0 bn,kak is given by 

(1 + xc{x'^))f{xc{x'^)) — ao 



x{l — xc{x'^)) 

Proof. This follows from the fundamental theorem of Riordan arrays since 



1 X 



1 + 1 + 



{c{x ),xc{x )) 



Example 8. The g.f. of the B-transform of the Fibonacci numbers -F„ is given by 

1 + xc(x^) xc(x^) 
x(l — xc(x^)) 1 — xc(x^) — x^c(x^)^ 

This follows since 



1 — X — x"^ 

and Fq = 0. This g.f. may be simplified to 

1 - + xa/1 - 1 - + x\/l - 4x2 

1 - 2x - 5x2 + iQ^3 - (1 - 2x)(l - 5x2) 

By solving the equation 

1 

u = 



x2 



1 + X — X2c('x2) 



□ 



□ 
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we see that this g.f. may be expressed (using Eq. (3)) as the continued fraction 



l-3x 



2x' 



l + x 



X 



1 - 



X 



X 



which shows that the the Hankel transform of the B-transform of the Fibonacci numbers is 
(-2)". 

We have defined the B matrix using the Riordan array {y^,^^)- This matrix is 
associated with the Chebyshev polynomials of the second kind Un{x) (it is the coefficient 



array of Un{x/2)). The matrix 
first kind T„. We have 

Proposition 9. We have 



l—x^ x_ 

l+a;2 ' 1+2 



is related to the Chebyshev polynomials of the 



X 



X 



1 + ' 1 + 



(1 + x,x) ■ {1 + X, x)*. 



Proof. We have 



L ■ [1 + x,xy = L ■ {1 + X, x)~^ ■ {l+x,x) ■ {1 + X, xY 
1 



L- 



,x 



l + x 

- X X 
2 



(1 + ■ (1 + x.xf 

1 



X^ 



,x 



(1 + x,x) ■ {1 + x,x) 



X 



X 



x^ 



x^ 



(1 + x,a;) ■ (1 + 



We can decompose (1 + x, x) ■ (1 + x, x)* as the sum of two matrices 



/ 1 

1 1 

110 

110 

1 1 

1 1 

\ \ \ \ \ \ \ 



\ 



/ 1 

110 

110 

1 1 

1 1 

1 

\\ \ \ \ \ \ 



\ 



(10) 



□ 



9 



which is the sum of (1 + x,x) and a shifted version of (1 + x,x). To obtain B we multiply 

-1 



by U 



L+a;2' l+2;2 



. This gives us, once again 





1 














. 


.. \ 




( 


1 











. 


.. \ 




1 


1 











. 









1 


1 








. 






2 


1 


1 








. 









2 


1 


1 





. 






3 


3 


1 


1 





. 




+ 





3 


3 


1 


1 


. 






6 


4 


4 


1 


1 


. 









6 


4 


4 


1 


1 . 






10 


10 


5 


5 


1 


1 . 









10 


10 


5 


5 


1 . 




v 


















v 















where the first member of the sum is the Riordan array 



X 



X 



1 + x'^ 1 + 



(1 + x,x) 



X 



X 



(1 +X)(1 + ,t2)' 1 + X 



L. 



Theorem 10. Let bn = Ylk=o bn,kCi-k where bn,k is the {n, k)-th element ofM. Then 



E 

k=0 



[0'n-2k + 0'n-2k+l) 



where we have extended an to negative n by setting a^n = o-r. 
Proof. We have seen that B has general term 



n 



Thus the B transform of a„ is given by 

( ( ["ziij 



n+l 

E 

A:=0 ^ 

which can also be written as 

n+l 



n 



n 



Qfe 



Qfe 



n 

LtJ 



n 

LtJ 



n 



a-k 



since O'^ ■ (^j) = O'^ ■ (|n^|). We can also write this as 



n+l 

E 

fc=0 



n 



(1-0 



n-fc+l^ 



n 

LIJ 



(11) 
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Now note that 

LfJ 



\ k/ ^^^'"^^ ~^ 0,n-2k+l) — ( ^ j {0'n-2k + 0'n-2k+l) + ( ^ j {0"n-2k + On-2fc+l) 

fc=0 ^ ^ fc=0 ^ ^ fc=LtJ+i 

LfJ 



fc=0 ^ ^ fc=LtJ+i 



By gathering similar terms in the above expression, and considering the cases of n even 
(n\2 = 0) and n odd, we arrive at 



(«n-2fe + an-2/c+l) = ^ (a„_2fc+l +2an-2fc + an-2fc-l) + [^^VS = 0] ^ ^ (ao + fli ) • 

(12) 

By considering the separate sums for k even and k odd in Eq. (11), extending to negative n 
and gathering terms we find that also 

["-^ J 

&n = ^ (a„_2fc+i + 2a„_2fc + a„_2fc-i) + [n\2 = 0] ^^^^^ (ao + fli). 



Thus we have the following equivalent expressions: 

n+l 

bn = bn^kO-k 

k=0 

n 



□ 



fc=0 

n-2A; + Q'n-2k~l 



fc=0 

fc=0 ^ 2 / ^ + ^ ^ j=0 i=l 



4 Symmetric Toeplitz plus Hankel matrices 

We now recall result Proposition 2.1 from [1], which we state in the language used above. 
Proposition 11. [1, Proposition 2.1]. Let {an)Z=-oo be a sequence with an = a-n and let 



k=0 



[0'n-2k + 0'n-2k+l)- 
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Also let H = ibi+j)i,j>o be the Hankel matrix of {bn)n>o o-nd A = {a-i^j + aj+j+i)jj>o be 
the Toeplitz plus Hankel matrix associated to (a„)"__oQ. Finally let L be the matrix with 
{n,k) — th term {^rJ^^)- Then 

H = L-A-L\ (13) 

An immediate consequence of this is that 

A = L-^H{L')-^ = L-'H{L-y. 

If now H has an LDU decomposition H = C - D ■ C'' then we obtain an LDU decomposition 
for the symmetric Toephtz plus Hankel matrix A: 

A = L-^ ■ C- D ■ C ■ {L-y, 



or 

A = (L-^C) ■ D ■ {L-^Cy. 
Example 12. We continue our example with the Fibonacci numbers. Thus let 

n+l 



bn,kFk 

k=0 



1 - + xVl - 



;i-2x)(l-5x2) 

For this sequence, we have the following LDU decomposition oi H = {bi+j)ij>o. 



/ 1 3 7 17 39 91 

3 7 17 39 91 207 

7 17 39 91 207 475 

17 39 91 207 475 1075 

39 91 207 475 1075 2445 

91 207 475 1075 2445 5515 

V ; ; : : ; ; 





1 














. 




/ 1 














. 






1 



















3 


1 











. 







-2 











. 






3 


1 
















7 


2 


1 








. 










-2 








. 






7 


2 


1 













12 


6 


2 


1 





. 













-2 





. 






12 


6 


2 


1 










39 


13 


7 


2 


1 


. 
















-2 


. 






39 


13 


7 


2 


1 







91 


33 


15 


8 


2 


1 . 



















-2 . 






91 


33 


15 


8 


2 


1 


v 
















v '■■ 














V 















(14) 



Here, the first matrix C of the product is the inverse of the coefficient array of the orthogonal 
polynomials for which the sequence bn is the moment sequence. These polynomials are 
specified by 



Pn{x) = xPn^lix) - P„_2(X), P^^x) = 1, Pi(x) 



X 



S,P2{X) 



X 



-2x-l. 
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We have 



( 1 

















2 


1 














3 


1 


1 











5 


2 


1 


1 








8 


3 


2 


1 


1 





13 


5 


3 


2 


1 


1 



and thus 



A 



















V 


f 


1 














. 


.. \ 


/ 1 




2 


1 











. 









3 


1 


1 








. 









5 


2 


1 


1 





. 









8 


3 


2 


1 


1 


. 









13 


5 


3 


2 


1 


1 . 







V 














•• / 


V 



1 


2 


3 


5 


8 


13 


2 


2 


4 


6 


10 


15 


3 


4 


5 


9 


14 


23 


5 


6 


9 


13 


22 


35 


8 


10 


14 


22 


34 


56 


13 


16 


23 


35 


56 


89 





















-2 



















-2 



















-2 



















-2 
























/ 1 



















2 


1 
















3 


1 


1 













5 


2 


1 


1 










8 


3 


2 


1 


1 







13 


5 


3 


2 


1 


1 



/ V 



We note that the matrix L~^C in this case is "almost" a Riordan array, in that it is the 
Fibonacci "sequence-array" ( ^_J_^2 , x) with general term [k < n]F„„fc+i, shifted once with 
a first column of Fn+2 pre-pended. 

Example 13. We take the example of the Jacobsthal numbers 

{-IT 



2" 

y 



X 



\x 



1 — X ~ 2x^ 



We note that this is the the element corresponding to r = 2 of the family of sequences with 
n-th term given by 



X 



\x 



1 — X — rx^ 



k=0 



n — k — 1 
k 



where the Fibonacci numbers correspond to r = 1. Thus we let 



n+l 



bn,kJk 



k=0 



13 



Then the g.f. for bn is given by 

(c(x), xc(x^)) 



2^^ , (l+^)(Tr^)-0\ _ v/r^4^^ + 3(l-2a;) 



x{l — x) 



This is equivalent to the expansion 



2(2 - 9x + 10x2) 



l-3a; 



X 



X 



X 



1 

from which we deduce that the Hankel transform of the B-transform of the Jacobsthal 
numbers is (—1)". Using Eq. (9), we can also write the g.f. of bn as 



1 — X X 



1 + ' 1 + 
and hence we have 



-1 /(1 + x) 



{l+x){l-2x) 



X 



1 — X X 



1 + x^ ' 1 + x^ 



1 



1 - 2x' 



I — n I — n ^ L 2 J / 



k=0 k=0 2 

The Hankel matrix H for bn has LDU decomposition LDL^ as follows: 



H 



( 1 


3 


8 


21 


54 


138 


3 


8 


21 


54 


138 


350 


8 


21 


54 


138 


350 


885 


21 


54 


138 


350 


885 


2230 


54 


138 


350 


885 


2230 


5610 


138 


350 


885 


2230 


5610 


1408J 





1 














. 


• ^ 


/ 1 














. 


■■ ^ 




1 














. 


.. \ 




3 


1 











. 







-1 











. 






3 


1 











. 






8 


3 


1 








. 










-1 








. 






8 


3 


1 








. 






21 


9 


3 


1 





. 













-1 





. 






21 


9 


3 


1 





. 






54 


24 


10 


3 


1 


. 
















-1 


. 






54 


24 


10 


3 


1 


. 






138 


64 


27 


3 


1 


1 . 



















-1 . 






138 


64 


27 


3 


1 


1 . 




v 
















v 














v 














• / 



In this case, the matrix £ is a Riordan array, equal to 
C = 



VI - 4x2 + 3(1 - 2x) 



2(2 - 9x + 10x2) 



xc(x 



1 - 3x + 2x2 



X 



x^ 



x" 



14 



Here, y i+t^^ ; T+ii^j coefficient array of the family of orthogonal polynomials given 

by 

P„(x) = xP„_i - P„_2(a^), Po{x) = 1, Pi{x) = X - 3, P2ix) = - 3x + 1. 

The B-transform of the Jacobsthal numbers J„ is thus the moment sequence for this family 
of orthogonal polynomials. Finally, we have 



1 — X X 



Vl -4x2 + 3(1 -2x) , 2, 
-,xc(x^) 



l + x2'l + xV V 2(2 -9x + 10x2) 
1 — X X \ /I — 3x + 2x2 X 



1 + X2' 1 + x2 

.X 



1 + X2 ' 1 + X^ 



1 - 2x 



Thus the Toeplitz plus Hankel matrix A associated to the Jacobsthal numbers J„ has LDU 
decomposition 



A 



/ 1 2 4 8 16 32 

2 3 6 12 24 48 

4 6 11 22 44 88 

8 12 22 43 86 172 

16 24 44 86 171 342 

32 48 88 172 342 683 
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